RESCEU-02/11 



f{R) Gravity and its Cosmological Implications 

Hayato Motohashi Alexei A. Starobinsky and Jun'ichi Yokoyama ^'^ 
^ Department of Physics, Graduate School of Science, 
^ _ The University of Tokyo, Tokyo 113-0033, Japan 

^ ' ^ Research Center for the Early Universe (RESCEU), 

, ^ . Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan 

^ L. D. Landau Institute for Theoretical Physics, Moscow 119334, Russia 
Q I ^ Institute for the Physics and Mathematics of the Universe(IPMU) , 

u: The University of Tokyo, Kashiwa, Chiha, 277-8568, Japan 

^\ Abstract 
O 

5-H ■ We have investigated the evolution of a homogeneous isotropic background of the Universe and 
c/2 . 

. inhomogeneous subhorizon matter density perturbations in viable f{R) models of present dark 
energy and cosmic acceleration analytically and numerically. It is found that viable f{R) models 

\^ ' generically exhibit recent crossing of the phantom boundary wde = — 1- Furthermore, it is shown 

^ ! 

■ that the growth index of perturbations depends both on time and wavenumber. This anomalous 

o _ 

growth may explain properties of the observational matter power spectrum from the SDSS data 
and can also partially counteract the spectrum suppression by massive neutrinos making larger 
values of the total sum of neutrino rest masses possible. 

X 
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I. INTRODUCTION 



Although the standard spatially flat A-Cold-Dark-Matter (ACDM) model can explain 
cosmic acceleration and is consistent with current observational data[l], the observed value 
of the cosmological constant term is much smaller than any other energy scale known in 
physics. On the other hand, we are sure that "primordial dark energy (DE)", which is 
responsible for inflation in the early universe [2-4], is not identical to the cosmological 
constant, in particular, it is not stable and eternal. Hence, it is natural to seek nonstationary 
models of the current DE, too. 

f{R) gravity is one of those dynamical DE models in which the Hilbert-Einstein action 
is modified and generalized by incorporating a new phenomenological function of the Ricci 
scalar R, f{R). This theory provides a self-consistent and nontrivial alternative to the 
ACDM model. It contains a new scalar degree of freedom dubbed "scalaron" in Ref. [2]. 
The existence of this additional degree of freedom imposes a number of conditions on viable 
functional forms of f{R) for R^ Rq and up to curvatures R in the center of neutron stars: 



where the prime denotes the derivative with respect to the argument R and Rq is a parameter 
of the order of the present Ricci curvature. Furthermore, f{R) should satisfy the stability 
conditions: 



Specific functional forms satisfying all these conditions have been proposed in Refs. [5-7]. 

In this paper, we have carried out numerical calculations of the evolution of both a 
background space-time and density fluctuations for the particular f{R) model introduced 
in Ref. [7]. As a result, we have found the phantom boundary crossing at an intermediate 
redshift z < 1 for the background space-time metric and an anomalous behaviour of the 
growth of density fluctuations. 
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FIG. 1: Evolution of the equation-of-state parameter wde{z) for the effective dark energy. 
II. BACKGROUND UNIVERSE 

We adopt the following action with a three-parameter family of f{R) models: 

^ = / d'x^f{R) + S^, (3) 

iOTrCr J 



f{R) = R + XRs 



(4) 



where n, A and Rg are model parameters and Sm is the action of the matter content. We 
can derive field equations from the action (3) and rewrite them as 



(5) 
(6) 



where F{R) = f'{R)- Working in the spatially flat Friedmann- Robertson- Walker (FRW) 
space-time with the scale factor a{t), we find 



3H^ = SnGp - 3{F - l)H^ + ]^{FR - f) - 3HF 
2H = -SnGp - 2{F - 1)H - F + HF, 



(7) 
(8) 



where H = a/a is the Hubble parameter and p is the energy density of matter. 

From Eq. (6), we can express the effective energy density and pressure of dark energy as 



^ttGpde = -3HRF' + 3{H^ + H){F-1)- ^(/ - R), 



87rG(pDE + Pde) = 2H{F -1)~HF + F, 



(9) 
(10) 
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TABLE I: -B(O) for various model parameters. 



n X 


B{0) 


2 0.95 


2 OQ X ID" 


-1 


2 4 




-4 


2 8 


6.07 X 10" 


-5 


3 0.73 


1.86 X 10" 


-1 


3 2 


1.34 X 10" 


-3 


3 3 


1.35 X 10" 


-4 


4 0.61 


1.73 X 10" 


-1 


4 1 


1.33 X 10" 


-2 



respectively where R = 12H'^ + 6H. We define the DE equation of state parameter wde by 
the ratio wbe = -Pde/pde, 



DE 



= -1 + 



2H{F -1)-HF + F 



(11) 



Pde -3HRF' + 3(H^ + H) (F - 1) - (/ - R) /2 

For the appropriate initial condition corresponding to the existence of cosmic inflation in 
the past, f — R acquires an asymptotically constant value f — R = —XRs at high redshift. 
In this regime the evolution of the Universe is the same as that obtained from the Einstein 
action with a cosmological constant A(oo) = XRs/2. 

The late-time asymptotic de Sitter solution has a curvature R — Ri which is obtained as 
the maximal solution of the equation 2/(i?i) = namely, 

"1 + (n + l)r^ 



a(r) =r + 2X 



= 0, 



(12) 



'I _|_ ^2^n+l 

where r = Ri/Rg. It is obvious that the Minkowski space-time, r = 0, is one of the solutions. 
The stability condition of this future de Sitter solution[8, 16], f'{Ri) > Rif"{Ri), imposes 
the following constraint on r, 

(1 + r2)[(l +r2)"+i - 2nXr] 



r > 0, 



(13) 



2nA[(2n + l)r2 - 1] 

which is stronger than any other constraint discussed above. For each n, we can find r, which 
marginally satisfies Eq. (13) and gives the minimal allowed value of A. Numerically we find 
{n, r^in, Ainin) = (2, 1.267, 0.9440), (3, 1.041, 0.7259), and (4, 0.9032, 0.6081) for n < 4. 
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We numerically solve the evolution equation (8) using Eq. (7) to check the numerical 
accuracy and taking ti as the moment of time when the matter density parameter was 
Qi — ISnGpi/ {IGnGpi + XRg) = 0.998. We determine the current epoch by the requirement 
that the value of Q takes the observed central value Qq = 0.27. Then Rg is fixed in such a 
way that the current Hubble parameter Hq = 72 km/s/Mpc is reproduced. 

Figure 1 depicts the evolution of wbe as a function of redshift z [14]. Phantom crossing 
is manifest there. As expected, wde approaches — 1 = constant as we increase A for fixed 
n. For A = Amin, deviations from wbe — —i are observed at ~ 5% level in both directions 
for z < 2 independently of n. Such behaviour of wde is well admitted by all the most 
recent observational data[l]. The average value of wj^-^ over the interval < 2; < 1 to which 
all BAO and most of the SN data refer is very close to —1. Moreover, in this range the 
behaviour of wde for A — Amin is well fitted by the CPL fit [9] wde(^) ~ wq + Waz/ {1 + z) 
with {n,Wo,Wa) = (2,-0.92,-0.23), (3,-0.94,-0.22), and (4,-0.96,-0.21), respectively. 
These values of wq and Wa he very close to the center of the 68% and 95% CL eUipses for all 
combined data in Fig. 13 of Ref. [1]. This phantom crossing is not peculiar to the specific 
choice of the function (4) but a generic one for models which satisfy the stability condition 
F' > 0. 

We have also calculated the quantity B{z) = {f / f){dR/dlnH) at the present time. 
The results are presented in Table I. 



III. DENSITY FLUCTUATIONS 



Armed with the evolution of the Universe background, we proceed to the investigation 
of the evolution of matter density fluctuations, 6, in f{R) gravity. In the subhorizon limit, 
the evolution equation is derived as [10, 11] 

S + 2HS - AnGospd = 0, (14) 

Equation (14) reduces to the correct evolution equation for all wavenumbers for the ACDM 
model in the Einstein gravity where F = 1. The time and /c-dependence of the effective 
gravitational constant Ggfr changes the evolution of density fluctuations. 
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FIG. 2: Evolutions of the growth index 7(2) and the effective gravitational constant Gcs{z) for 
n = 2 and A = 1. 

First, let us consider the gravitational growth index, 7(2;), which is an important quantity 
helping to distinguish different modified theories of gravity. It is defined through 



din (5 



or 7(^) 



(16) 



c? In a log f2m 

In the standard ACDM model, it takes a practically constant value 7 = 0.55. However, 
it evolves with time in modified gravity theories in general. We also note that 7(2;) has a 
nontrivial ^-dependence in f{R) gravity, since density fluctuations with different wavenum- 
bers evolve differently. Therefore, this quantity is a useful diagnostic to distinguish DE in 
modified gravity from the ACDM model in the Einstein gravity. We present the evolution 
of 7(z) together with that of Gcs/G for different values of k in Fig. 2. 7(2;) takes a constant 
value identical to the ACDM model in the early high-redshift regime because f{R) gravity 
is indistinguishable from the Einstein gravity plus a positive cosmological constant then. It 
gradually decreases with time, reaches a minimum, and then increases again towards the 
present epoch. Current constraints on the growth index[12] are not strong enough to detect 
any deviation from the ACDM model and/or to obtain new bounds on f{R) DE models, 
but future observations may reveal its time and wavenumber dependences. 

Second, we focus on the additional late-time transfer function for linear matter perturba- 
tions arising in f{R) gravity. The evolution equation (14) can be solved analytically in the 
high- curvature regime when the scale factor evolves as a(t) oc t^^^ and F takes the asymptotic 
form F ~ l-2nA (i?//?^)"^""^ = 1 - (i?/i?c)"^"\ where N = 2n and Rc = i?,(2nA)i/(2n+i)_ 
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FIG. 3: Additional transfer function C{k) in f{R) gravity. 



The two independent solutions of Eq. (14) in this regime read 



Skit) = 5ife ( — 



X 2F1 



±5 



33 ±5 



;1± 



4(3iV + 4)' 4(3A^ + 4)' 2(3A^ + 4) 



-;-3 



{N + i)e ft 



aim 



2N+8/3"' 



;i7) 



in terms of the hypergeometric function[13]. Hereafter, we consider the upper sign solution 
only, because the other solution corresponds to the decaying mode and is singular at t — )■ 0. 
Then, the solution behaves as 

-1+733 

(18) 



-j and 5fc(t) — ^ 5i,,C{k) 



respectively. The transfer function, C{k), is given by 



C{k) 



r 1 + 



2(3Af+4) 



33 \ 



2(3Af+4) J 



-r f-\ _i_ _5W33_\ -p f _5W33_^ 
J- 4(37v+4) J ^ \^4{3Af+4) J 



3{N+l)k' fSRctf 



am. 



2 \ N+2 



-5+%/33 
4(3iV+4) 



(19) 



where U = sinh~^ a/(1 — We have confirmed the additional transfer 

function C{k) numerically as Fig. 3. 

Finally, we proceed to compare the results of our numerical calculations for matter power 
spectrum with observational data from SDSS DR7[17]. The left panel of Fig. 4 shows the 
fitting of the matter power spectrum for the ACDM model. The vertical axis is the ratio of 
the power spectrum normalized by the power spectrum in ACDM model with bias parameter 
b = 1.25. The data with error bars are LRG samples in SDSS DR7. As far as we use a 
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FIG. 4: Power spectrum with SDSS data. Left: Power spectrum in the ACDM model normalized 
by -Pacdm(& = 1.25). Right: Power spectrum in f{R) gravity for n = 2 and A = 1 with total 
neutrino mass up to 0.3eV normalized by Pacdm(^ = l-l)- 

constant bias, we see that observational data are increasing with wavenumber compared with 
the theoretical expectation line in the ACDM model. Furthermore, if we include massive 
neutrinos, the fitting becomes worse because their free streaming suppress the matter power 
spectrum. We have found that f{R) gravity enhances the matter power spectrum and 
produces a better fit to the SDSS data than the ACDM model without introducing a scale 
dependent bias or a nonlinear bias, see the right panel of Fig. 4. There is also a possibility 
to admit a larger total sum of neutrino restmasses compared to the ACDM model since the 
anomalous growth of perturbations partially counteracts their suppression by free streaming 
of massive neutrinos [15]. 



IV. CONCLUSION 



In this paper, we have numerically calculated the evolution of both a homogeneous 
isotropic background and matter density fluctuations in a viable f{R) DE model based 
on the specific functional form proposed in Ref. [7]. We have found that viable f{R) gravity 
models of the present DE and accelerated expansion of the Universe generically exhibit phan- 
tom behaviour during the matter-dominated stage with crossing of the phantom boundary 
W£)E = — 1 at redshifts z < 1. More exactly, this behaviour is characteristic for all f{R) DE 
models that have f"{R) > and a stable future de Sitter epoch and which approach the 
Einstein gravity sufficiently fast for R ^ Rq, under the condition that the gravitational con- 
stant G in the Einsteinian representation of the field equations (4) is normalized to its value 
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measured in laboratory experiments (i.e., for it! ^ Rq, too). The predicted time evolution of 
wde has qualitatively the same behaviour as that has recently obtained from observational 
data[18]. However, it is important that the condition of stability, or even metastability, of 
the future de Sitter epoch strongly restricts the possible deviation of wde from —1 by several 
percents in these models. Thus, the DE phantomness should be small, if it exists at all, that 
agrees well with the present observational data. 

As for the density fluctuations, we have also investigated the growth index 7(/c, z) of 
density fluctuations and have presented an explanation of its anomalous evolution in terms of 
the time dependence of Geff- Note that this evolution is characteristic for all f{R) models in 
which the scalar particle (scalaron) becomes relativistic (A:^ > rnl{R)a^) at recent rcdshifts. 
Since 7 has a characteristic time and wavenumbcr dependence, future detailed observations 
may yield useful information on the validity oi f{R) gravity through this quantity, although 
current constraints have been obtained assuming that it is constant both in time and in 
wavenumber [12] . 

We have also numerically confirmed a shift in the power spectrum index for larger 
wavenumbers which exceed the scalaron mass during the matter-dominated epoch[13], while 
for smaller wavenumbers, fluctuations have the same amplitude as in the ACDM model. 
Once more, the future de Sitter epoch stability condition bounds a possible increase in den- 
sity fluctuations for cluster scales (compared with the ACDM model) by ~ 40% for n > 2. 
This enhancement in matter power spectrum can explain the observational data from SDSS 
and allows for some increase in the total sum of neutrino restmasses, as compared to the 
standard ACDM model. 
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